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ABSTRACT 


A  model  has  been  developed  describing  the  expected  costs 
of  delivery  and  delay  per  period  when  demand  for  a  certain 
repair  part  comes  from  more  than  one  production  line  at  a 
NARF.  This  model  extends  the  earlier  work  by  McMasters  and 
Davidson  on  a  scheduled  delivery  model.  The  objective  of 
the  model  development  was  the  determination  of  the  optimal 
number  of  periods  between  deliveries.  This  is  that  number 
N  which  minimizes  the  expected  costs  of  delivery  and  delay 
per  period.  Unfortunately,  no  simple  closed  form  expression 
for  optimal  N  as  a  function  of  the  other  parameters  could  be 
obtained.  As  a  consequence,  parametric  analyses  of  the  cost 
function  were  conducted  to  determine  optimal  N  and  its 
behavior  when  other  parameters  were  varied. 
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I .  INTRODUCTION 


A.  BACKGROUND 

In  October,  1979,  Naval  Air  Rework  Facility  (NARF) 

Alameda  was  the  first  such  facility  to  have  its  wholesale 
supply  support  provided  by  a  Naval  Supply  Center  (NSC) 
Oakland.  This  was  the  result  of  a  recommendation  from  the 
DoD  Material  Distribution  study  that  wholesale  supply  sup¬ 
port  of  NARF's  near  the  NSC's  at  Oakland,  San  Diego,  and 
Norfolk  be  provided  by  the  NSC's  instead  of  by  the  Naval 
Air  Stations  (Ref.  1]. 

The  question  of  providing  supply  support  for  the  local 
NARF,  with  no  degradation  of  that  support,  became  a  very 
real  concern  to  the  NSC's.  The  centralization  of  materials 
at  the  supply  centers  has  increased  the  distances  to  these 
wholesale  customers;  therby  escalating  the  costs  associated 
with  moving  material.  Not  only  does  the  transportation  cost 
increase,  but  more  importantly,  delays  are  created  by  ex¬ 
tended  travel  time  and  longer  picking  times.  These  delays 
cost  money  since  components  may  wait  longer  for  parts. 
Maintenance  time  must  be  used  to  gather  components  awaiting 
parts  and  store  them  together  until  required  parts  are  re¬ 
ceived.  Because  of  this,  valuable  production  resources  may 
be  diverted. 

Grant  [Ref.  2]  attempted  to  quantify  the  production  delay 
costs  when  repair  parts  are  not  immediately  available  when 
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required.  He  was  unable  to  establish  a  firm  relationship 
between  delivery  times  and  delay  cost  due  to  the  lack  of 
historical  data  at  the  NARF.  He  did  determine  that  the 
slowest  delivery  sets  the  pace  of  the  repair  action  and 
should  be  used  to  determine  the  production  delay  costs. 

Under  the  assumption  that  such  delay  costs  will  even¬ 
tually  be  quantified,  McMasters  [Ref.  3]  developed  three 
direct  delivery  models  as  a  first  step  in  determining  the 
best  way  to  support  the  NARF.  Each  alternative  considered 
involved  a  single  production  (overhaul)  line  placing  a  single 
demand,  at  any  instant  of  time,  on  the  supply  system.  Each 
demand  was  for  a  repair  part  required  by  a  component  being 
reworked.  The  demand  for  this  part  was  considered  as  a 
Bernoulli  trial  with  a  fixed  known  probability  of  demand  (p) 
for  each  component  inducted. 

Due  to  the  analytic  complexities  of  McMasters'  models  no 
closed  form  optimization  formulas  were  possible.  Therefore, 
the  expected  total  costs  were  used  in  calculating  the  optimal 
delivery  schedules.  Davidson  [Ref.  4]  performed  a  parametric 
analysis  on  these  three  models  concentrating  primarily  on  the 
probability  of  demand  and  production  line  delay  costs.  She 
discovered  that  there  was  little  difference  in  the  optimal 
expected  costs  for  each  model.  The  ^jst  practical  alterna¬ 
tive  appeared  to  be  to  make  sche  i  deliveries. 

Minimization  of  the  expected  costs,  where  the  total  cost 
was  the  sum  of  transportation  costs  and  delay  costs,  was  the 
goal  of  the  work  done  by  McMasters  and  Davidson.  The 


transportation  cost  was  defined  to  be  a  fixed  charge  per 
delivery.  Delay  cost  was  a  fixed  charge  per  component  per 
unit  of  time  delayed  after  it  was  identified  as  needed  for 
a  given  component's  overhaul. 

B.  PURPOSE 

A  key  assumption  of  the  McMasters -Davidson  work  was  that 
only  one  production  line  created  the  demand  for  a  given 
repair  part.  Since,  in  reality,  more  than  one  line  might 
need  a  given  part,  this  assumption  should  be  relaxed  and  the 
effects  of  that  relaxation  investigated.  Therefore,  this 
thesis  will  extend  their  scheduled  delivery  model  to  the 
consideration  of  multiple  production  lines. 

C.  PREVIEW 

Chapter  II  presents  the  assumptions  and  the  derivation  of 
the  formulas  describing  the  expected  costs  per  period  for 
any  number  of  production  lines  needing  a  given  part.  Chap¬ 
ter  III  then  follows  the  procedure  of  Davidson  and  determines 
the  optimal  delivery  schedule  for  a  variety  of  cost  and 
probability  parameter  values  and  numbers  of  production  lines 
needing  the  part.  Chapter  IV  will  summarize  the  findings  of 
the  thesis  and  present  conclusions  and  recommendations  for 
further  areas  to  be  investigated. 
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A.  INTRODUCTION 

A  NARF  inducts  components  for  repair  and,  as  the  first 
step  of  the  process,  identifies  those  parts  which  need  to  be 
replaced.  Requisitions  are  then  prepared  and  transmitted  to 
the  supply  system.  That  system  picks  the  item  from  inventory 
and  delivers  it  to  the  NARF.  Upon  receipt  of  the  required 
parts  from  the  supply  system,  the  repair  actions  are  com¬ 
pleted  and  the  item  is  returned  to  the  system,  either  for 
immediate  installation  or  as  a  shelf  spare. 

All  of  the  earlier  studies  considered  alternative  trans¬ 
portation  systems  for  delivering  a  given  requisitioned 
repair  part  from  the  central  stock  point  to  the  customer. 

They  also  attempted  to  minimize  the  sum  of  expected  trans¬ 
portation  costs  and  expected  customer  delay  costs.  McMasters 
[Ref.  3]  considered  two  basic  direct  delivery  scheduling  schemes, 
scheduled  delivery  and  unscheduled  delivery.  Under  scheduled 
delivery,  a  deliver/  truck  would  be  dispatched  to  the  NARF 
every  N  periods  of  time  and  would  deliver  all  requisitions 
that  had  been  submitted  up  to  that  time  by  the  NARF.  Under 
unscheduled  delivery  the  truck  would  make  a  delivery  as  soon 
as  K  requisitions  had  been  received.  Davidson  found  that 
the  minimum  expected  total  costs  of  delivery  and  delay  were 
essentially  the  same  for  either  scheme.  This  then  allows  a 
selection  of  the  scheme  which  is  the  most  desirable  from  other 


aspects  of  the  problem.  Scheduled  delivery  is  considered  to 
be  more  desirable  since  it  is  easier  for  the  NSC  to  plan 
and  provides  the  NAP.F  with  a  known  delivery  time.  Scheduled 
delivery  also  represents  the  way  many  stock  points  presently 
operate  their  local  delivery  systems. 

B.  ASSUMPTION'S  MADE  IN  GENERALIZING  FROM  PREVIOUS  MODELS 

In  generalizing  to  a  multiple  production  line  system, 
some  assumptions  from  the  earlier  model  must  be  modified. 
McMasters'  model  referred  to  a  time  period  as  "the  time 
between  component  inductions  cn  the  production  line."  For 
the  multiple  production  line  environment,  different  produc¬ 
tion  lines  may  have  different  periods  between  inductions. 
However,  within  a  group  of  similar  lines,  such  as  aircraft 
engines,  the  induction  quantities  for  a  giver,  quarter  have 
been  observed  to  be  fairly  close  in  value.  In  addition, 
since  such  schedules  are  negotiated  ahead  of  time  for  a  given 
quarter  at  workload  planning  conferences,  it  is  possible 
to  obtain  identical  schedules  for  several  lines.  Thus,  this 
thesis  will  assume  equal  length  time  periods  for  ail  produc¬ 
tion  lines. 

As  in  the  case  of  the  earlier  model,  transportation  costs 
will  be  considered  as  a  fixed  charge  per  ro- -c  ip  from  the 
supply  center  to  the  receiving  facility.  is  not  an 

unreasonable  assumption  since  all  parts  for  similar  lines 


are  usually  received  at  a  central  point  at  the  NARF  and  then 
dispersed  internally  to  the  various  production  lines. 

Ke  will  also  retain  the  assumption  that  Chambers  [Ref.  5. 
made  in  his  model.  It  was  that  all  requirements  were  homo¬ 
geneous;  that  is,  no  consideration  was  given  as  to  priority 
of  individual  production  lines  over  other  product icr.  lines. 

Delay  costs  are  assumed  to  be  different  for  each  produc¬ 
tion  line.  We  will  consider  only  time  dependent  delays.  The 
fixed  costs  associated  with  a  delay,  which  includes  the  cost 
of  removing  the  component  from  the  production  line  and  gather¬ 
ing  associated  parts,  will  not  be  considered. 

The  probability  of  demand  for  a  given  repair  part  by  a 
given  line  is  assumed  to  be  known  in  advance.  However,  these 
probability  values  will  be  assumed  to  be  different  for  each 
line. 

C.  DETERMINISTIC  DEMAND 

If  a  demand  from  a  customer  occurs  once  every  time  period 
with  certainty,  it  is  said  to  be  a  deterministic  demand.  Let 
CT  be  the  cost  of  one  round  trip  from  a  supply  center  to 
the  NARF.  If  a  truck  is  dispatched  every  time  a  demand  is 
received  and  processed,  the  cost  of  delivery  for  each  unit 
is  C,p.  If,  on  the  other  hand,  the  truck  waits  until  k  units 
have  been  demanded  and  processed,  the  average  delivery  cost 
per  unit  is 


CT/k  • 


(2.1) 
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This  is  also  the  average  delivery  cost  per  period.  If  the 
truck  waits  until  it  is  full  (suppose  it  has  a  capacity  of 
n  units)  then  the  delivery  cost  per  period  is  minimized  at 


CT/n  •  (2.2) 

While  the  k  units  are  accumulating,  they  are  creating 
delay  costs  for  the  NARF.  If  the  truck  waits  for  k  units  to 
be  accumulated  and  the  delay  cost  for  one  unit  for  one  time 
period  is  CD,  the  total  average  delay  cost  per  period  is 

CD(k-l)/2  *  (2.3) 

To  verify  Equation  (2.3),  assume  one  unit  is  needed  every 
time  period.  If  the  truck  waits  for  k  units  to  accumulate 
it  will  not  leave  until  (k-1)  time  periods  after  the  first 
demand.  During  this  tii.ie  the  units  ordered  but  not  delivered 
have  caused  delay.  Specifically,  the  first  unit,  ordered 
in  the  first  period,  will  be  delayed  (k-1)  periods;  the  second 
unit,  ordered  in  the  second  time  period,  will  be  delayed 
(k-2)  periods;  etc.  Only  the  kth  unit  ordered  in  period  k 
will  have  no  delay.  The  total  waiting  time  in  periods  then 
is 

(k-1)  +  (k-2)  +  (k-3)  +  ...  +  1  +  0  ,  (2.4) 

which  can  be  written  as  k(k-l)/2.  When  this  is  multiplied 
by  the  delay  cost  per  period,  the  result  is  the  total  delay 


k 
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cost 


The  average  delay  cost  per  period  is  obtained  by  dividing 
by  the  number  of  units,  k. 

CD(k-l)/2  (2.6) 

By  adding  the  average  shipping  cost  and  delay  costs  per 
period,  the  total  average  cost  for  k  units  becomes 

TAC(k)  =  CT/k  +  CD(k-l)/2  .  (2.7) 

TAC(k)  is  a  discrete  function  since  k  must  be  integer  valued. 
As  a  consequence,  minimization  requires  that  finite  differ¬ 
ences  be  used.  The  optimum  k  is  that  value  such  that 

TAC(k-l)  >  TAC(k)  <  TAC ( k+1)  (2.8) 

or  equivalently  the  largest  k  such  that 

TAC(k)  -  TAC(k-l)  <  0  ;  (2.9) 

or  the  smallest  k  such  that 

TAC(k)  -  TAC (k+1)  <  0  .  (2.10) 

Using  Equation  (2.7),  inequality  (2.9)  becomes 
CT/k  +  CD(k-l)/2  -  (C^k-1)  +CD(k-2)/2]  <  0  ,  (2.11) 
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which  reduces  to 


k(k-l)  v  2cm/'cD  •  (2.12) 

This  final  relationship  allows  a  very  simple  iterative 
computation  step  to  be  made  starting  with  k  =  1  and  repeated 
until  the  largest  k  is  found  which  still  satisfies  Equation 
(2.12)  . 

D.  DERIVING  THE  EXPECTED  COST  FUNCTION  FOR  MULTIPLE 

PRODUCTION  LINES 

McMasters  [Ref.  3]  derived  a  formula  for  the  expected 
delivery  and  delay  costs  based  on  a  single  demand  source 
over  an  infinite  time  horizon.  In  his  derivation,  he  repre¬ 
sented  the  round  trip  cost  of  a  delivery  by  CT,  as  was  done 
for  the  model  above.  That  notation  will  be  retained  in  this 
thesis.  We  will  assume  that  will  be  incurred  each  time  a 

delivery  takes  place.  However,  if  a  demand  does  not  occur 
during  an  interval  of  N  periods,  there  will  be  no  cost 
incurred  due  to  a  delivery  cancellation. 

The  delay  costs  were  defined  above  as  being  incurred  by 
the  NARF  as  a  result  of  not  having  a  needed  part  available 
at  the  time  it  is  required.  Several  elements  could  be  con¬ 
sidered  as  part  of  «-hat  cost.  These  elements  include  the 
cost  associated  with  putting  the  component  aside  (such  as 
documentation,  parts  ordered  and  putting  it  on  the  shelf) , 
labor  costs  due  to  work  stoppage,  inventory  costs  and  cost 
associated  with  part  non-availability  to  the  customer.  For 


the  purpose  of  this  thesis  the  delay  cost  will  be  denoted 


by  CD. 

The  expected  costs  per  period  for  a  single  prccuct.cn 
line  as  derived  by  McMasters  (Ref.  3]  is: 


N  infl  M  CT  PC.  (  N-  1) 

ecp(n)  =  (i-d-p)N)  t-i!lLL_LL_ELJ.]  [_I  *—5, - 

d-p)N  N  1 


(2.13) 


A  comparable  form  can  be  derived  for  the  multiple  production 
line  case.  The  steps  of  the  derivation  follow  those  of 
Reference  3. 

First  the  expected  total  delays  for  a  given  production 
line  must  be  determined.  The  delays  are  a  function  of  the 
number  of  different  configurations  that  demands  can  take 
over  the  N  periods  between  deliveries.  As  an  example, 
suppose  N  =  2,  then  there  are  four  possible  configurations. 
First,  a  demand  occurs  in  the  first  period  and  none  occurs 
in  the  second,  resulting  in  a  delay  of  one  period.  Next,  a 
demand  does  not  occur  in  the  first  period,  but  does  occur 
in  the  second.  This  results  in  a  delay  of  zero.  The  third 
configuration  for  N  =  2  is  a  demand  in  both  the  first  and 
second  periods .  A  total  delay  of  one  period  results  from 
the  first  period’s  demand  (zero  delay  for  the  second  period's 
demand)  resulting  in  a  total  delay  of  one  period. 

The  total  number  of  configurations  for  production  1,  when 
exactly  x^  demands  occur  is 
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Therefore,  the  total  number  of  configurations  which  can 
occur  over  N  periods  and  have  at  least  one  demand  is 


n 


l 


N  N  N, 

l  nx  =  I  (/> 

x„=l  xl  x.  =  l  xl 


N,-l 

r. 


(2.15) 


The  expected  total  delay  time  for  the  1th  production 
line  associated  with  the  N.  configurations  can  be  determined 
in  the  following  manner.  First  consider  only  those  configura¬ 
tions  having  exactly  x^  demands  where  x^  >_  1 .  The  probability 
of  each  configuration  is: 


N„  -x . 


P<VSt' 


■  Pi 


(2.16) 


The  number  of  configurations  having  a  demand  in  period 
1  <  iZ  <  is 

N  -1 

m,  =  (x^)  •  (2.17) 

It  is  important  to  realize  that  m^  is  independent  of  i,. 

Those  demands  occurring  in  period  i^  will  have  to  wait 
until  N^-i^  periods  have  passed  before  the  parts  are  delivered. 
The  total  of  all  delays  for  these  configurations  having  x^ 


demands  is  therefore: 


T 


TD,  (X, ,N  ,) 

•  V,  i 


NZ  N  -1 

•»  X  ,  1  *•  A. 

1^1  X 


N  -1  Ni 

■  J.  <W 

4  1=1 

X 


N. ( N  „  - 1 ) 

X  X 

- 2 — 


N ,  -1 

X 


! 


(2.18) 


From  Equations  (2.16)  and  (2.18)  it  follows  that  the  total 
expected  delays  for  the  1th  production  line  over  all  values 
of  x.  are: 

i 


N 


ETD^(N) 


l  TD. (x, ,N,)P(x. :N.  ) 


x^l 


N 


l  N1(N,-1)  N^-l  xz 


N ,  -X . 

V.  K 


*iml 


N.  (N-l)p,  Ni  N.-l  >1,-1  N  ,-x . 

Xx  v,  r  /  x  /i  _  \  4  *• 


L(*-i)pi  ^-p4> 


x;=l  ^ 


(2.19) 


Notice  that  the  summation 


N, 


Xt*l 


N 


l 

Z 


-1 


(1-P,) 


(2.20) 


is  the  sum  of  the  binomial  probabilities  over  all  values  of 
x^  _>  1  and  its  value  is  1.  Because  of  this,  Equation  (2.19) 
reduces  to 

N  (N  -Dp 

ETD£(N£)  =  — - ^ - -  •  (2.21) 

In  developing  the  expected  costs  per  period  over  all 
production  lines,  consideration  must  be  given  to  the  fact 
that  a  delivery  will  occur  at  the  end  of  N  periods  if  there 
has  been  at  least  one  demand  during  that  time  interval.  In 
the  case  of  multiple  production  lines  this  also  means  that 
at  least  one  line  had  a  demand  during  N.  Using  two  produc¬ 
tion  lines  to  begin  the  development,  there  are  three  mutually 
exclusive  events  that  can  occur  in  N  periods  which  would  place 
a  demand  on  the  system.  The  first  is  that  production  line 
one  has  at  least  one  requirement  and  production  line  two  has 
no  requirements.  The  expected  cost  per  period  for  the  N 
periods  would  be 

C  PiCD  <N-1) 

ETC (N)  =  [  1- ( l-px)  N )  (l-p2)N(1I  + - ±- - )  .  (2.22) 

V 

Here  [l-(l-p1)  ]  represents  the  probability  that  production 

N 

line  one  has  at  least  one  demand  and  (l-p2)  represents  the 

probability  that  production  line  two  has  no  demand  in  the  N 

periods.  Another  event  is  when  production  line  one  has  no 

N 

requirements  (with  probability  (1-p^)  )  and  production  line 

N 

two  has  at  least  one  demand  (with  probability  [ 1— ( 1 — p2 )  ])• 
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The  expected  cost  would  be: 


C  ^2<"D 

ETC (N)  =  [l-(l-p2)N]  (l-p1)N[1I  + - L - ]  .  (2.23) 

Finally  production  line  one  has  at  least  one  requirement 
and  production  line  two  has  at  least  one  requirement.  The 
expected  cost  would  be: 

p,C  (N-l)  p  C  (N-l) 

ETC (N)  =  [l-(l-Pl)N]  [l-(l-p2)N]  [-£  + - 1 -  + - - j 

(2.24) 

The  total  expected  costs  from  the  two  production  lines  for 
the  N  periods  is  the  sum  of  the  costs  associated  with  the 
three  possible  events. 

C  P1CD  (N-1} 

ETCT(N)  =  u-(i-Pl)N]  (i-p2)N[1I  + - L- - j 

C  P2^D 

+  [1-(1-P2)N](1-Pi)N[^  + - \ - ] 

N  N  CT  PlCD1(li*1) 

+  [i-(i-Plrj  ii-d-p2)N)  i-i  +  — f - 

P2C  <N-1) 

+  — h — 1 

(l-(l-Pl)N(l-p2)N]CT  U-(l-pL)N]p1CD^(N-l) 

=  -  +  - 


T 


If  three  production  lines  were  under  consideration  the 
total  expected  costs  for  N  periods  would  be  as  follows: 


N 


N 


N  i 


[l-(l-p.)"(l-p,)"(l-p,)”]CL.  :i"(1"pl)N]plCD 
ETCt(M)  =  - ± - j— - - - -  +  - 2 - - - 


N 


[l-(l-p2)Hp2CD  ( N-l ) 


tl- (l-p3)  LN]p3CD^(N-l) 


2.26) 


For  L  production  lines,  the  total  expected  delay  cost  function 
for  the  first  N  periods  would  be: 


[l-(l-Pl)N(l-p?)N--*(l-pT)N]C 
ETCt(N)  =  - - - jp - - - - 

[l-(l-Pl)N]PlCD  ( N-l) 

*  - 2 - 1 - 

[l-(l-p2)N]p2CD  (N-l) 

+  2 

[1-(1-Pl)N)PlC  (S-l) 

+  ...  +  - 2 - - -  .  (2.27) 


The  generalized  form  for  a  finite  number  of  production  lines, 
L,  for  N  periods  is 


ETCT(N) 


[1  - 


L 

n 

i=i 


L  U-(l-p/]PeCD  (N-l) 

l  - 2 - 1 - 

i»l 


(2.28) 
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Now  there  are  actually  four  mutually  exclusive  events 
that  can  take  place  if  there  are  two  production  lines  under 
consideration.  The  fourth  event  is  that  no  demands  will 
occur  in  either  production  lines  one  or  two.  The  probability 
of  this  event  is  ( 1-p^) N ( l-p2) N .  Because  of  this  event,  a 
renewal  theory  approach  which  considers  a  sequence  of  many 
intervals  of  N  periods  is  appropriate.  A  renewal  occurs  if 
a  delivery  occurs.  The  first  three  events  of  the  two  produc¬ 
tion  line  model  correspond  to  a  renewal  in  the  first  N  periods 
whereas  the  fourth  corresponds  to  renewals  in  later  N  period 
intervals.  In  this  latter  case  consideration  must  be  given 
to  what  the  expected  total  costs  will  be  in  the  second  set 
of  N  periods.  For  the  two  production  lines,  the  form  of  the 
expected  total  costs  over  the  second  N  periods  depends  on 
four  mutually  exclusive  events  that  can  take  place.  The 
first  is  that  neither  production  line  one  nor  two  places  a 
demand  in  the  first  N  periods  but  production  line  one  has  a 
demand  during  the  second  N  periods.  The  expected  costs  per 
period  are: 

C  ^l^D 

m  m  N  N  x 

ETC  ( 2N)  =  (l-p1)W(l-p2)N[l-(l-p1)N]  (l-p2)N[^+ - ±4 - ]. 

(2.29) 

The  second  event  is  the  same  as  the  first  except  that  line 
two  has  a  demand  and  line  one  does  not.  The  expected  costs 
per  period  are: 
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ETC ( 2N) 


■  “p  i  j 


p  C  (N-l) 

r  T  ,  L2 

2N  ^  4 


]  • 


The  third  event  is  when  both  lines  one  and  two  only  have 
demands  during  the  second  N  periods.  The  expect.ec  delay  cc: 
per  period  are 


ETC ( 2N)  =  {l-p1)N(l-p2)N[l-(l-p1)N}  ll-(l-p2)N,l 

r  p,C_  (N-l)  p-Cn  (N-l) 
i  u.  *  e>2 

•4  +  — \ - +  — h — i  • 


(2.31) 


The  fourth  event  is  that  no  demands  occur  during  the  second 
N  periods.  Expected  costs  per  period  when  a  renewal  spans 
2N  periods  is  the  sum  of  the  above  three  events. 


C  P1CD  (N'1) 

etct(2N)  =  (i-p1)Na-p2)N[l-(l-p1)Nni-P2)NlI|  + - 3 - 3 


f  ■,  I  N  /  i  .  N  r  i  ^  »  N  *  »  ■%  *  N  r  T 

+  (1-P1)  (1-P2)  [  1-  ( 1-P2)  ]  (1-P1)  [55?  + 


P-)C_  (N-l) 
2 

4 


+  ( 1-Pl ) N ( 1-P2) N [1-{1-pl)N] [1-(1-p2)N] 
c  p^  (N-l)  P2CD  (N-l) 

+ - T - + - T - i 


This  simplifies  into  Equation  (2.32) 
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ETC_ ( 2N) 


( l-p,) N(l-p2) N  Ct!1-( 1-Pi)N(1-p2)NJ 

-  3  -  -  N  1 

p,  C„  (N-l)  [l-(l-t>.  j'!] 

1  W1 

.  - - - 

P2CD  (N-l)  U-(l-p2)N] 

+  - 1 - ^ - .  .  (2.22) 

If  we  apply  the  generalized  form  for  more  than  two  production 
lines,  the  expected  total  costs  per  period  associated  with 
the  second  set  of  N  periods  is 


ETC*  ( 2N) 


t  (l-o  )X  CTil*  a”p..’ 
•  1  25  1 


,=  1 


N 


p,Cn  (N-l)  ll-(l-pf )-] 
1 - 


(2.23) 


If  no  demands  occur  during  the  first  and  second  sets  of 

N  periods,  the  third  set  of  N  periods  must  be  considered. 

Assuming  there  has  been  no  demand  in  the  first  and  second 

intervals,  then  the  probability  associated  with  this  occur- 

2N  2N 

rence  would  be  < 1— p ^ )  *  < 1 — p 2 )  ‘  in  the  two  production  lines 
case.  The  expected  total  costs  per  period  due  to  delay  when 
3N  periods  occur  before  a  renewal  is 


ETCtv3N) 


2m  2M  m  m 

(i-p-, )  d-p0?  cT::-(i-p1)J,'i-p2)3'] 

_ ; _ . _ •  — : - 


N 


PoC  (N-l)  [1-d-p,) 

£.  U  *- 


V.  A.  »  ^  Tt  I 


The  generalized  form  for  more  than  two  production  lines 
for  the  third  interval  becomes 


ETCt(3N) 


^  2N  ^ 

:  (i-p  )  cmu-  :  (I-?  )  : 

=1  x  -.=i 

3  N 


L  ptCD  (N-l)  [1- ( 1-p  )N] 
—  *. 

■  =1  2 


Generalizing  to  k  intervals  of  N'  periods,  the  total  expected 
costs  per  period  for  L  production  lines  is: 


ETCT(kN) 


-  n  0  v  (k-l)N  f1  - 
..  (1-p  /  C_U-  .. 

■>1  '  ■  1  .:  =  1 


( 1-P  >  ] 
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\  _ L 
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■1  =  1 


(N-l)  [1-d-p  J 
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(2.36) 
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Finally,  the  total  expected  costs  per  period  over  all 
future  periods  can  be  obtained  by  summing  over  all  possible 
k  values: 


ETCt(N) 


7  (1-p.)  {k_i)K  c  [1-  7  (1-p 
1  =  1  '•  ■  1  =  1 


k=l 


k 


N 


+ 


p  C D  (N-l)  [1-  (1-p  .  )  K] 
2 


(2.37) 


The  summation  term  in  Equation  (2.37)  can  be  rewritten  as 
follows : 


,  "  ‘1-P,) 
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(k-1)  N 


k  =  l 
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1=1 
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and  since 
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Therefore  Equation  (2.37)  reduces  to 


ETCt(N) 


n  d-p.)!:] 

;.=i 


(i-P,) 


N 


L 

U-  n  (i-p, > c, 

-  ,  4 

4  =  1 


N 


L  P(CD  (N-l)[l-[l-p/’] 

+  7  - - - 5 - ;■  .  (2.38) 

4  =  1  i 


E.  DETERMINATION  OF  OPTIMAL  N 

Since  N  can  take  on  only  integer  values,  it  is  necessary 
to  use  f.inite  differences  to  determine  the  optimal  value  for 
N.  Optimal  N  is  that  value  of  N  which  satisfies  the 
following: 


ETCt(N-1)  >  ETCt(N)  <  ETCt(N+1)  .  (2.39) 

This  is  the  same  as  saying  that  optimal  N  is  the  largest 
value  of  N  such  that 

1ETCt(N)  =  ETCt(N)  -  ETCt(N-1)  <  0  .  (2.40) 

However,  if  one  attempts  to  evaluate  the  form  of  AETC(N),  it 
becomes  quickly  apparent  that  the  result  is  more  complex  than 
ETC,p  ( N )  itself.  Therefore,  the  determination  of  optimal  N 
in  the  next  chapter  consists  of  evaluating  ETC^,  (N)  for 
increasing  N  as  long  as  inequality  (2.40)  is  satisfied. 


F.  COST  FUNCTION  IF  PROBABILITY  OF  DEMAND  EQUALS  ONE 


For  purposes  of  comparison  in  the  next  chapter  the 
special  case  of  Equation  (2.38)  when  P]^  =  P2  =  •••  =  P^  =  1 
will  be  useful.  In  that  case.  Equation  (2.38)  reduces  to: 


ETCt(N) 


L 

V  x  - 


(2.41) 


which  is  a  simple  extension  of  Equation  (2.7)  and  a  delivery 
will  occur  at  the  end  of  the  first  N  periods. 
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III.  OPTIMIZATION  ANALYSIS 


A.  RELEVANT  VARIABLES  IN  DETERMINING  AN  OPTIMAL  N 
Determination  of  an  optimal  N  is  centered  around  the 

influence  of  the  cost  of  transportation  from  the  supply 
point  to  the  NARF  and  return,  the  cost  due  to  delay,  the 
probability  of  demand  and  the  number  of  production  lines. 

The  intent  of  the  optimization  analysis  will  be  to  show  how 

N  varies  with  respect  to  p„,  C  ,  and  1.  The  transportation 

4  U£ 

cost  CT  will  be  fixed  at  $100  throughout  the  analysis  since 
it  is  independent  of  the  number  of  production  lines. 

B.  GENERAL  BEHAVIOR  OF  THE  COST  CURVES 

Figures  3.1  through  3.5  are  graphical  representations  of 
the  expected  total  costs  per  period  and  its  two  component 
parts,  transportation  costs  and  delay  costs  per  period. 
Because  the  domain  of  the  cost  function  is  integer  valued, 
each  point  on  a  graph  represents  the  cost  at  that  N  value. 
The  points  have  been  connected  for  better  visualization  of  a 
cost  function's  behavior.  Figures  3.1  through  3.3  are  for  a 
single  production  line  and  show  the  effect  of  increasing  C^. 
As  the  value  of  CD  is  increased,  the  value  for  the  optimal  N 
decreases.  In  all  cases  the  value  for  the  CD  curve  starts 
at  zero  when  N  =  1  and  the  slope  of  the  curve,  as  the  value 
of  CD  increases,  becomes  steeper.  The  ETCT(N)  curve 
approaches  the  curve  asymptotically. 


30 


Figures  3.1,  3.4  and  3.5  illustrate  the  impact  from  the 
number  of  production  lines.  Both  the  transportation  and 
delay  cost  components  are  affected.  In  these  figures  the 
pi  and  CD  values  are  set  equal;  p£  =  0.1  and  CD  =  $10. 

As  the  number  of  lines  increases  the  cost  term  increases. 
It  is  most  significant  for  N  =  1.  The  influence  from  the 
number  of  lines  then  decreases  rapidly  as  N  increases. 

The  Cp  cost  term  is  always  zero  when  N  =  1  regardless  of 
the  number  of  production  lines.  However,  as  N  increases  it 
increases  at  an  increasing  rate.  When  there  are  more  than 
one  production  lines,  the  rate  of  increase  is  faster  because 
all  of  the  lines  experience  delays.  The  additive  effect  is 
however  less  than  linear  with  the  number  of  lines. 

Figures  3.6  through  3.8  present  the  effects  of  various 
parameter  changes  on  ETCT  for  a  two  production  lines  case. 

In  these  figures  Cn  and  p.  are  fixed  at  $10  and  0.1,  respec- 

l>i  l 

tively.  Figure  3.6  then  sets  N  =  2  and  shows  the  effects 

that  changing  and  p^  have  on  ETC^,(2)  .  Figures  3.7  and 

3.8  show  similar  results  for  N  =  3  and  N  =  5,  respectively. 

For  a  given  Cn  value  the  total  costs  appear  to  increase 

with  increasing  N.  However  it  is  not  true  for  all  p^  values 

as  Figures  3.9  and  3.10  show.  These  figures  have  fixed 

L2 

at  $10  and  $20,  respectively. 

Figures  such  as  3.9  and  3.10  can  be  used  to  determine 
optimal  N  as  a  function  of  p^ •  Optimal  N  corresponds  to  the 
minimum  ETCT  curve  for  any  given  P2  value.  For  example, 
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Figure  3.6  Expected  Total  Cost  Curves  as  a  Function  of 
P2  for  Two  Production  Lines  When  CD  is 

Varied  and  C_  =  $10,  p,  =0.1  and  $  =  2 


Figure  3.7  Expected  Total  Cost  Curves  as  a  Function  of 
P2  for  Two  Production  Lines  When  CQ  is 

Varied  and  C_  =  $10,  p.  =0.1  and  =  3 


Figure  3.9  shows  that  N  =  2  is  never  opimal  for  the  values 
of  the  parameters  used  in  our  study,  that  N  =  4  is  optimal  for 
p-,  from  approximately  0.9  3  to  1.0,  that  N  =  5  is  optimal  for 
from  0.58  to  0.93  and  N  =  6  is  optimal  from  p^  =  0  to 
0.58  (actually  N  =  6  is  not  optimal  over  the  whole  range 
because  larger  N  values  yield  lower  cost  turves) .  Figure  3.11 
shows  the  plot  of  optimal  N  as  a  function  of  for  the  data 
from  Figure  3.9  plus  additional  curves  for  N  =  7  through  12. 

As  is  expected,  the  optimal  value  of  N  decreases  with  in¬ 
creasing  P2 ,  the  probability  of  a  demand  in  any  period  from 
a  second  production  line. 

Figure  3.12  shows  the  plot  of  optimal  N  provided  by 

Figure  3.10  plus  additional  curves  for  N  =  7,  8,  and  9. 

Comparison  of  Figures  3.11  and  3.12  shows  the  impact  of 

changing  C  from  $10  to  $20.  As  expected,  ootimal  N  ae- 
2 

creases  with  increasing  C  . 

u2 

Figures  3.13  and  3.14  show  a  dramatic  decrease  in  optimal 

N  because  Cn  and  Cn  have  been  increased  to  $50. 

U1  u2 


Figure  3.15  shows  optimal  N  when  p.  =  0.5,  C_  =  $10 

i  u1 

and  Cn  =  $20.  This  figure  when  compared  with  Figure  3.14 
2 

shows  the  expected  result  that  optimal  N  decreases  as  p 
increases . 

Finally,  Figures  3.16  and  3.17  present  results  for  three 
production  lines  when  p1  =  p2  =  0.1  and  CD  =  $10  for  all 
three  production  lines.  The  effect  of  adding  a  line  can  be 
seen  by  comparing  Figure  3.17  with  3.11.  The  result  is  that 
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Figure  3.12  Optimal  N  for  Two  Production  Lines  When 
Cn  =  $10,  Cn  =  $20  and  p,  =0.1 
D1  d2 
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Figure  3.17  Optimal  N  for  Three  Production  Lines  When 
C„  =  =  $10  and  p.  =  p„  =  0.1 


Figure  3.17  looks  like  Figure  3.11  moved  to  the  left. 
Interestingly,  the  shift  in  break  points  is  approximate! 
0.1,  which  is  the  demand  probability  of  the  added  line. 
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IV.  SUMMARY  AND  CONCLUSIONS 

A.  SUMMARY 

A  model  has  been  developed  describing  the  expected  costs 
of  delivery  and  delay  per  period  when  demand  for  a  certain 
repair  part  comes  from  more  than  one  production  line  at  a 
NARF .  This  model  extends  the  earlier  work  by  McMasters  and 
Davidson  on  a  scheduled  delivery  model.  Several  assumptions 
were  made  to  facilitate  the  development.  First,  the  time 
periods  between  component  inductions  for  each  production 
line  were  considered  equal.  Next  the  transportation  costs 
were  fixed.  Thirdly,  the  time  dependent  delay  costs  were 
allowed  to  be  different  for  each  production  line.  The  proba¬ 
bilities  of  demand  for  each  production  line  were  also  allowed 
to  be  different. 

The  objective  of  the  model  development  was  the  determina¬ 
tion  of  the  optimal  number  of  periods  between  deliveries. 

This  is  that  number  N  which  minimizes  the  expected  costs  of 
delivery  and  delay  per  period. 

Unfortunately,  no  simple  closed  form  expression  for  optimal 
N  as  a  function  of  the  other  parameters  could  be  obtained. 

As  a  consequence,  parametric  analyses  of  the  cost  function 
were  conducted  to  determine  optimal  N  and  its  behavior 
when  other  parameters  were  varied. 
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B.  CONCLUSIONS 

The  optimization  analysis  showed  that  as  the  vaiu- 
Cq  was  increased  lor  eacn  production  tint-#  wLil*_'  j 

c 

the  number  of  lines  1  and  p  constant,  the  value  for  optical 
N  decreased.  When  the  value  of  CQ  is  held  constant  for  cr.e 
production  line  and  increased  for  another  under  the  above 
conditions  the  optimal  N  values  also  decreased.  The  decrease 
for  the  latter  case  is  not  as  severe  as  when  all  C^  values 
are  increased. 

When  the  probability  of  demand  is  increased,  holding  L 

and  CD  constant,  optimal  N  decreases.  For  small  probabiii- 

i. 

ties  of  demand  it  is  less  likely  that  repair  items  will  be 
required,  thereby  extending  the  number  of  periods  between 
deliveries.  However,  as  the  probabilities  increase  it  is 
more  likely  that  a  demand  will  be  made  in  the  earlier  periods. 

As  the  number  of  production  lines  increase,  the  optimal 
value  of  N  also  decreases.  Here  again  the  decrease  corres¬ 
ponds  to  the  fact  that  the  likelihood  of  a  demand  from  the 
system  has  increased. 

This  study  assumed  all  demands  were  treated  equally  by 
the  NARF.  In  reality,  this  is  not  the  case.  Each  demand 
submitted  has  a  priority  attached  to  it;  however,  within  a 
given  production  line,  a  repair  part  usually  has  the  same 
priority  each  time  it  is  requested.  Thus,  to  illustrate  the 
impact  of  this  priority,  it  is  only  necessary  to  change  a 
production  line's  CD  value  to  a  higher  value  to  reflect  its 
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increased  importance .  At  present  the  priority  system  is  ar. 
implicit  way  to  distinguish  between  the  fact  that  different 
lines  have  different  CD  values  because  the  h'ARF  cannot  com¬ 
pute  actual  values. 

This  thesis  was  only  able  to  explore  the  parametric 
influences  to  a  limited  extent.  Additional  analyses  sugges 
themselves.  For  example,  why  dc  the  break  points  in  Figure 
3.17  appear  to  have  shifted  by  approximately  0.1  from  those 
in  Figure  3.11  when  one  production  line  was  added.  The 
shift  of  0.1  was  the  same  as  the  probability  of  demand  for 
the  added  production  line. 

Perhaps  more  important  is  the  need  to  evaluate  the  rang 
of  parameter  values  fcr  which  N  =  1  since  that  corresponds 
a  delivery  as  soon  as  a  demand  occurs.  This  compares  to  th 
current  twice  daily  delivery  service  provided  by  NSC  Gaklar. 
to  NARF  Alameda.  The  results  cf  such  an  analysis  could  be 
used  tc  determine  if  current  actual  parameter  values  were 
comparable  or  net.  If  not,  then  perhaps  less  frequent 
service  is  worth  considering. 
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